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Abstract — In this paper we propose and discuss coherent terahertz 
sources based on charge density wave (plasmon) amplification in two 
dimensional graphene. The coupling of the plasmons to interband 
electron-hole transitions in population inverted graphene layers can lead 
to plasmon amplification through stimulated emission. Plasmon gain 
values in graphene can be very large due to the small group velocity 
of the plasmons and the strong confinement of the plasmon field in the 
vicinity of the graphene layer. We present a transmission line model for 
plasmon propagation in graphene that includes plasmon dissipation and 
plasmon interband gain due to stimulated emission. Using this model, we 
discuss design for terahertz plasmon oscillators and derive the threshold 
condition for oscillation taking into account internal losses and also losses 
due to external coupling. The large gain values available at terahertz 
frequencies in graphene can lead to integrated oscillators that have 
dimensions in the 1-10 ^m range. 
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I. Introduction 

Graphene is a single two dimensional (2D) atomic layer of carbon 
atoms forming a dense honeycomb crystal lattice yj. The electronic 
properties of graphene have generated tremendous interest in both 
experimental and theoretical arenas \T\, f4l, fS], f6], fT|, fSl, f9l, 
[12]. The energy dispersion relation of electrons and holes with zero 
(or close to zero) bandgap results in novel behavior of both single- 
particle and collective excitations (l], |2|, 14^, JSJ, ||6J, |7|, [81, |91, 
1121 . The high mobility of electrons in graphene has prompted theo- 
retical and experimental investigations into graphene based ultra high 
speed electronic devices such as field-effect transistors, pn-junction 
diodes, and terahertz devices (7), fM, El, ID, 0, Qll, fT3l , 
1141. Negative conductance at terahertz frequencies under population 
inversion conditions has been predicted by Ryzhii et. al. in 1131 , 
n4 :l and the potential for graphene based terahertz oscillators was 
also suggested there. It has been shown that the frequencies of the 
charge density waves (or plasmons) in graphene at moderate carrier 
densities ( 10^ - 10^^ cm"^) are in the terahertz frequency range |4|, 
l5l . 1121 . Plasmons in graphene are strongly coupled to the interband 
electronic transitions. The zero bandgap of graphene leads to strong 
damping of the plasmons at finite temperatures as plasmons can 
decay by exciting interband electron-hole pairs (2), (3), 1121 . Plas- 
mon emission with the accompanying decay of electron-hole pairs 
has also been experimentally observed in graphene |15|. Recently, 
the author showed that plasmons in graphene can experience very 
large gain values (exceeding 10* cm^^) at frequencies in the 1- 
10 terahertz range under moderate population inversion conditions 
due to the stimulated emission of plasmons [12|. This process is 
depicted in Fig. \T\ Plasmon gain in graphene is similar to the optical 
gain resulting from the stimulated emission of photons in III-V 
semiconductor interband lasers at optical and infrared frequencies. 
The main difference is that the slow group velocity of plasmons in 
graphene at terahertz frequencies and the strong confinement of the 
electromagnetic field associated with the plasmons near the graphene 
layer results in plasmon gain values that are much larger than the 
typical gain values in semiconductor interband lasers. 

In this paper, we propose and discuss terahertz oscillators based 
on plasmon amplification in graphene. We present a transmission line 
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Fig. 1. (LEFT) Energy bands of graphene showing stimulated absorption of 
plasmons. (RIGHT) Population inverted graphene bands showing stimulated 
emission of plasmons. 



model for plasmon propagation in graphene that includes plasmon 
dissipation due to intraband scattering and plasmon interband gain 
due to stimulated emission. Using this model, we discuss terahertz 
plasmon gain in graphene, present designs for practical terahertz 
oscillators, and derive the threshold condition for oscillation taking 
into account both intrinsic losses and losses due to coupling power 
out of the device. A critical parameter that emerges from this 
analysis is the characteristic impedance of the plasmon transmission 
line. The threshold condition is shown to depend on the ratio of 
the external impedance and the characteristic impedance of the 
plasmon transmission line. The transmission line model presented 
here is not just relevant for plasmon oscillators but can also be 
used to model the high frequency electrical behavior of graphene 
based devices and interconnects. We use this model to obtain the 
impedance of a graphene strip under population inversion conditions. 
The electromagnetic energy of plasmons is confined within very small 
distances of the graphene layer in the vertical direction. In Section lVTl 
we discuss techniques to confine plasmons in the lateral direction as 
well using sub-micron scale dielectric waveguiding structures in order 
to minimize propagation losses. The large gain values predicted for 
plasmons in graphene can enable very compact integrated terahertz 
amplifiers and oscillators. 

Plasmons in metals and semiconductors from THz to optical 
frequencies have been the subject of much attention in the last few 
years 1,1 6J , 117J . Plasmon losses in metals are generally large in the 
visible and near-IR frequency range 1161 . 1171 . The strong confine- 
ment of the plasmon field in the vicinity of the metal surface near 
the plasmon resonance frequency makes amplification of plasmons 
difficult to achieve with a gain medium placed outside the metal. 
Plasmons in graphene are particularly interesting and unique since 
the material that supports the plasmons can also provide plasmon 
gain. The strong confinement of the plasmon field in the vicinity of 
the graphene layer then becomes an advantage and not an obstacle 
for achieving plasmon amplification. 

II. Plasmons in Graphene: The Random Phase 
Approximation (RPA) 

In this Section, we summarize one of the methods that has been 
used previously to describe plasmons in graphene and is based 
on the random phase approximation (2), (3), 1121 . In the next 
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Section, we will present a transmission line model for plasmons in 
graphene that is more relevant for device applications. In graphene, 
the valence and conduction bands resulting from the mixing of the 
Pz-orbitals are degenerate at the inequivalent K and K' points of 
the Brillouin zone 1 1 1. Near these points, the conduction and valence 
band dispersion relations can be written as fTI . 
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where v is the velocity of the electrons and holes and equals 10 
cm/s JT]. We assume different Fermi levels for conduction and 
valence electrons to allow for non-equilibrium electron-hole popu- 
lations. The electron and hole densities are given by the expressions. 
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The factor of 4 in the front accounts for spin degeneracy and the two 
valleys at K and K' \ l\, and /(£c(k) - Efc) and /(£„(k) - Ef^) 
are the Fermi distribution functions of the conduction and valence 
electrons with Fermi energies Ef^ and -E/„, respectively. The 
wavevector k is measured from the K(K') point. The complex 
propagation vector q((jj) of the plasmon mode of frequency uj is 
given by the equation, e(q, i^j) — 0, where e{q,uj) is the longitudinal 
dielectric function of graphene (2), (3), 1121 . In this plasmon mode, 
the charge density (not the number density) associated with the 
electrons and the charge density associated with the holes oscillate 
in-phase. In the random phase approximation (RPA) e(q,6t;) can be 
written as L18J, 



e(q,<.^) = l-y(q)n(q,c^) 



(4) 



Here, ^(q) is the bare 2D Coulomb interaction and equals e^/2ecx)g. 
Coo is the average of the dielectric constant of the media on either 
side of the graphene layer. n(q,a;) is the electron-hole propagator 
including both intraband and interband processes. For small wavevec- 
tors (i.e. for v\q\ < cj) the intraband and interband contributions to 
the propagator can be approximated as follows liT2l . 
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Here, q = |q|. In Equation the intraband contribution to the 
propagator is written in the plasmon-pole approximation that satisfies 
the f-sum rule 1181 . This approximation is not valid for large values 
of the wavevector q when uj{q) —> vq. However, in this paper we 
will be concerned with small values of the wavevector for which 
the plasmons have net gain. Plasmon energy loss due to intraband 
scattering has been included with a scattering time r in the number- 
conserving relaxation-time approximation which assumes that as a 
result of scattering the carrier distribution relaxes to the local equi- 
librium distribution 1 19|. The real part of the interband contribution to 
the propagator modifies the effective dielectric constant and leads to 
a modification in the plasmon frequency under population inversion 
conditions. The imaginary part of the interband contribution to the 
propagator incorporates plasmon loss or gain due to stimulated inter- 
band transitions. The real and imaginary parts of the propagator in 
Equations {5) and {6} satisfy the Kramers-Kronig relations. Equations 
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Fig. 2. A graphene layer connected with metal ohmic contacts on each end 
(top view). Note that the width W is not the physical width of the graphene 
layer - but the width of the region in which the plasmons propagate. 



(O and ^ can be used with Equation lO to calculate the real and 
imaginary parts of the plasmon frequency uj{q) as a function of the 
wavevector q. However, from the point of view of device design, it is 
more useful to assume that the frequency lj is real and the propagation 
vector q{uj), written as a function of lj, is complex. Equations 
and can be solved for the complex wavevector q. Since the charge 
density wave corresponding to plasmons has the form g'^ "'-*'*''^ jjjg 
imaginary part of the propagation vector corresponds to net gain or 
loss. 



III. A Transmission Line Model for Plasmons in 
Graphene 

A. Transmission Line Equations 

The procedure for obtaining the plasmon dispersion in graphene 
outlined in the previous Section does not explicitly account for the 
connection of the device with the outside world. In this paper, we 
present a transmission line model for the plasmons in graphene that is 
more suitable for device design applications. We assume a graphene 
layer of effective width W and length L, connected on each side 
with a metal ohmic contact to the outside world, as shown in Fig. |2] 
As discussed later in Section [Vll the width W is not necessarily the 
physical width of the graphene layer but the effective width of the 
region in which the plasmons can be confined. We also assume in 
this Section that population inversion can be achieved in the graphene 
layer. Details on how to achieve population inversion are discussed 
in Section IVIII We assume a plasmon wave of frequency uj and 
wavevector q propagating in the +2-direction for which variations in 
the charge density Q (units: Coulombs/m), current Iq (units: Amps), 
and the electrostatic potential V (units: Volts) in the plane of the 
graphene layer can be written as. 



Q{z,t) = -Q{q,oj) 



,ijj}e ^ + c.c. 



V{z,t) = -Viq,uj)e' 



+ c.c. 



+ c.c. 
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Outside the graphene layer the potential decays exponentially as 
g-i\y\^ where \y\ is the distance from the graphene layer. Charge 
conservation implies. 
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B. Electrostatic Capacitance and Quantum Kinetic Inductance 

The potential V{q,cj) can be related to the charge density using 
the electrostatic capacitance per unit length Ces{<1,^), 



V{q,Lu) 



where Cssiq,'^) i^, 



(11) 
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The expression above for the capacitance assumes a wide (Wq >> 1) 
graphene layer with no top or bottom metal plane. It also assumes 
that qL » 1. The condition qL » 1 will be assumed to be true 
throughout this paper It follows from Equations dlOl ) and i ll It that 
the current can be related to the potential as follows, 



ujCEs{q,i^) 



V{q,iu) 



(13) 



A variation in the potential generates a current that can also be 
expressed using the inductance per unit length L{q, cj) and the 
resistance per unit length R{q,u}) of the plasmon transmission line 
(as is done in standard transmission line analysis 125]), 



'iqV{q,uj) 



)L[q, 



R{q,Lo) 
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The inductance L{q, u)) is the sum of the kinetic inductance Lk {q, i^) 
and the magnetostatic inductance Lms{<1,^) 1201 , and the resistance 
R{q, uj) includes plasmon dissipation due to intraband momentum 
relaxation scattering of electrons and holes as well as plasmon gain 
due to interband stimulated emission. For all geometries considered 
in this paper, the magnetostatic inductance is found to be much 
smaller than the kinetic inductance, and will therefore be ignored. 
This has been shown to be true for carbon nanotubes as well |20|. 
The expressions for LK{q,i^) and R{q, uj) can be obtained by 
looking at the microscopic dynamic response of graphene to potential 
perturbations and this response is described by the electron-hole 
propagator n( 17,0;) [18], 



n{q,Lo)V{q,Lo) 



(15) 



The microscopic expressions for LK{q,(j->) and R{q,u>) can be 
obtained by comparing the right hand sides of Equations il4\ 
and (il5l. Unlike in the case of carbon nanotubes |20|, 121], the 
expression for the kinetic inductance in graphene is frequency and 
wavevector dependent. The dispersion relation for the plasmons can 
be obtained from the standard transmission line formula that relates 
the propagation vector and the frequency |25|, 

iLjCEsiq,i^) {LK{q,i^) + LMs{q,Lo)) + R{q,ui)] = q^ 

(16) 

Assuming LMs{q,LiJ ~ 0, Equation l ll6t simplifies to, 



CEs{q,i^) 



(17) 



The above equation for the plasmon dispersion is almost identical 
to the one obtained by setting e{q, uj) equal to zero in the random 
phase approximation (RPA). The only difference is that Coulomb 
interaction is expressed in a more general way in terms of the 
electrostatic capacitance in Equation l ll7t . Equation J17t . which is 
obtained using the transmission line model of the plasmons, thus 
correctly reproduces the plasmon dispersion relation. Given the 
plasmon dispersion relation, one may drop the explicit functional 
dependence on the wavevector q in all the expressions. 



C. Power Flow 

The power P(a;) carried by the plasmon wave can be found using 
the Poynting vector. 



dxdy [E{uj) X H*{u;)] 



dy V{u;) 



dt 



(18) 



(19) 



The second term in the parenthesis in Equation l ll8t is due to 
the displacement current and arises because unlike in a standard 
transmission line the plasmon field is not completely transverse 
to the direction of propagation. The displacement current term's 
contribution to the power is exactly one-half of the first term and has 
the opposite sign. Consequently, the expression for the total power in 
Equation il9\ is different by a factor of two compared to the standard 
transmission line formula [25|. In order to conserve power when the 
device is connected to an external circuit, it is convenient to define 
the transmission line current I{u>) as. 



(20) 



The current is the maximum current that can be generated by the 
plasmon wave in an external circuit consistent with the conservation 
of power This will be shown explicitly in Section |V] Using this 
definition of the current, the power carried by the plasmon wave can 
be written in the standard form. 



PH = -K{\/Hr(a;)} 



(21) 



D. Plasmon Transmission Line Impedance 

Finally, the expression for the characteristic impedance Zo{ijj) of 
the plasmon transmission line is defined as, 

V{uj) 



=2 



Lk(i^) + i R{lo)/lo 



Ces{lo) 



(22) 



TV. Terahertz Plasmon Gain in Graphene 

In this Section, we discuss the plasmon gain in graphene at 
terahertz frequencies. Few conditions must be met in order for the 
plasmons to have net gain. The condition K{g(a')} < lj/v must 
be satisfied by the plasmon dispersion relation in order to avoid 
direct intraband absorption of the plasmons by electrons and holes. 
The condition Efc — Efv > huj must be satisfied in order for the 
plasmons to have gain via stimulated emission at the frequency lo. In 
addition, for net gain the loss due to intraband (momentum relaxing) 
scattering of electrons and holes must be less than the gain due 
to stimulated interband transitions. This implies that the imaginary 
part of the propagator Il{q,Lj) must be positive (and the resistance 
R{q,u!) must be negative). We define the net plasmon energy gain 
g{u!) as — 2 3{g(aj)}. The net plasmon gain is the difference of the 
interband plasmon gain and the intraband plasmon losses. Figs. I3I6I 
show the plasmon dispersion and the net plasmon gain for different 
electron-hole densities at T=77K and T=300K. Plasmon dispersion 
(but not plasmon gain) in graphene has been previously reported in 
(21, O, |4|, |5|. The values used for v and too in our calculations 
are 10* cm/s and 4 to, respectively (assuming Silicon-dioxide on both 
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Fig. 3. Calculated plasmon dispersion relation in graphene at 77K is plotted 
for different electron-hole densities (n = p = 1, 2, 3, 4 X 10^" cm~^). The 
condition t<j(q) > vq is satisfied for frequencies that have net gain in the 
terahertz range. The assumed values of v and r are 10** cm/s and 0.5 ps, 
respectively. 



Fig. 5. Calculated plasmon dispersion relation in graphene at 300K is plotted 



for different electron-hole densities (n : 



P - 



1,1.5,2,2.5 X 10^ 



The condition t<j(q) > vq is satisfied for frequencies that have net gain in 
the terahertz range. The assumed values of v and r are 10* cm/s and 0.5 ps, 
respectively. 




Fig. 4. Net plasmon gain (interband gain minus intraband loss) in graphene at 
77K is plotted for different electron-hole densities (n = p = 1, 2, 3, 4 X 10^" 
cm~^). The assumed values of v and r are 10* cm/s and 0.5 ps, respectively. 



Fig. 6. Net plasmon gain (interband gain minus intraband loss) in graphene at 
300K is plotted for different electron-hole densities (n = p = 1, 1.5, 2, 2.5 X 
10^^ cm~^). The assumed values of v and t are 10* cm/s and 0.5 ps, 
respectively. 



sides of the graphene layer). The calculations assumed an intraband 
scattering time of 0.5 ps which is less than the experimentally 
observed scattering time of around 0.6 ps in mobility measurements at 
T=58K 1 12 1, |9|. In Ref. |12| it has been shown that the plasmon gain 
remains large even for scattering times as small as 0.1 ps. The electron 
and hole densities were assumed to be equal (i.e. Efc = —Ef^). It 
can be seen that the net gain approaches 4 x lO** cm^^ for carrier 
densities in the lO^^-S x 10^^ cm~^ range for frequencies in the 1-10 
THz range at T=77K and at T=300K. At low frequencies (less than 1 
THz), the losses from intraband scattering tend to dominate. At higher 
frequencies, the net gain becomes positive for electron-hole densities 
larger than a minimum threshold value. The calculated plasmon gain 
values in graphene are between 10 to 100 times larger than the typical 
material gain values of most Ill-V semiconductors in the optical and 
infrared frequency range [22J. This large difference is due to the 
small group velocity of plasmons in graphene at terahertz frequencies 
and is similar to the gain enhancement observed in photonic crystals 
near the photonic band edges due to the small group velocity |23). 
The calculated group velocity of plasmons in graphene is shown in 



Fig. |7] The large values for the plasmon gain could enable terahertz 
amplifiers and oscillators that are only few microns in length. Design 
of terahertz oscillators and the conditions necessary for achieving 
terahertz oscillation are discussed in the next Section. 

V. Terahertz Plasmon Oscillators: Design and 
Oscillation Condition 

We consider the terahertz plasmon oscillator depicted in Fig. [S] 
where a population inverted graphene strip of width W and length 
L is connected to an external circuit. The condition for achieving 
oscillation in a plasmon oscillator must take into account the energy 
coupled out of the device since this energy must also be compensated 
by the plasmon gain in addition to the intrinsic losses that arise due 
to plasmon loss. We assumed that the external impedance connected 
to the device (as seen from the two terminals of the device) is 
Zext{i^)- Zext{i^) also iucludcs the resistances associated with the 
ohmic contacts to the graphene layer. The external impedance could 
also represent an antenna structure for coupling the radiation out of 
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Fig. 7. The group velocity of the plasmons in graphene. normalized to the 
speed of light in free space, is plotted for different electron-hole densities 
(n = p = 1, 2, 3, 4 X lO^" cm-2) for T=77K. The assumed value of v is 
10* cm/s. 
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Fig. 8. (TOP) A schematic of a graphene terahertz oscillator coupled to an 
external load. (BOTTOM) Equivalent circuit model. 
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Fig. 9. The characteristic impedance Zoiyj) of a 0.5 /jm wide plasmon 
transmission line is plotted as a function of the frequency. 




Fig. 10. The threshold gain gth{^) needed for oscillation is plotted as a 
function of the cavity length L for different values of the parameter |r|. 



the device. We write the current and the potential associated with the 
plasmons in graphene as a superposition of forward and backward 
propagating waves using the complex time-harmonic transmission 
line notation 1251, 



V(z) = V+{uj)e''^ + V-{u;)< 



(23) 
(24) 



The reflection coefficient r(tj) for the plasmons at the two ends of 
the graphene strip is defined as, 

V-{lo) 



V+ico) 



(25) 



The boundary conditions imposed by the external circuit are the 
following, 

I{z = 0) = I{z = -L) 
V{z = 0) - V{z = -L) 



I[z = 0) 
These boundary conditions give. 



(26) 
(27) 



Zext{(^) + 2 Zo{uj) 



and 



i q L 



(28) 
(29) 



Equation ( \29\ gives the condition necessary for achieving terahertz 
oscillation. Since q = 3?{g} — j<?/2, the threshold gain required to 
achieve oscillation as well as the frequencies (or the wavevectors) 
of the oscillating plasmon modes can be obtained from Equation 
l |29t . The threshold gain gthii-^) needed to achieve oscillation can be 
expressed as, 

3t4<^) = T logT^^T^ (30) 



1 1 1 



The oscillation condition in Equation l l29t can also be obtained us- 
ing a slightly different approach. The frequency dependent impedance 
Zs{uj) of the graphene device looking in from its two terminals at 
z — and z — L can also be obtained using the transmission line 
inodel and is found to be. 



Zs{lo) = 2Zo{L0)r 



The oscillation condition in Equation l |29l l can then be written in the 
more familiar form as Zaiio) = ~Z^xt(y)). 
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Finally, it needs to be shown that the expression for current in 
Equation ( |2Ql l and the boundary conditions in Equations J26ll-J27t 
are consistent with the conservation of power between the oscillator 
and the external circuit. The net power leaving the plasmon oscillator 
from the two ends at z = and z = —L can be found using the 
expression in Equation l l21t . 

\V+{u)\^ 



2 X 



2Zo{u:) 

The net power delivered to the load is 

1/(^ = 0)'^ 



(i-irnr) 



(32) 



(33) 



Using Equation l |28K the above two expressions can be shown to be 
equal. 

A. Discussion 

The ratio of the characteristic impedance of the plasmon transmis- 
sion line to the external impedance determines the threshold gain. 
If the external impedance is much smaller or much greater than the 
impedance of the transmission line then the threshold gain required 
for oscillation is close to zero (since F ~ ±1). This means that 
oscillation can be achieved if the gain from stimulated emission 
overcomes only the internal transmission line losses but no power is 
coupled out of the device and the device has therefore little practical 
utility. If in some frequency range Zextiijj) ~ 2Zo{q,u->), then 
r(a') ~ 0, and the threshold gain is infinite. In this case, plasmons 
are not well confined in the cavity and oscillation cannot be achieved 
for this frequency range. Both situations - when the impedance of 
the plasmon transmission line is completely mismatched or perfectly 
matched to the external impedance - need to be avoided. A large 
capacitive external impedance can also be harmful as it can short 
the device at high frequencies and result in poor power output 
coupling efficiencies. Achieving oscillation as well as decent power 
output coupling efficiencies both therefore depend on the ratio of the 
external impedance and the characteristic impedance of the plasmon 
transmission line. The impedance of the plasmon transmission line 
can range from a few hundred Ohms to a few Kilo-Ohms. Fig. |9] 
shows the characteristic impedance Zoiuj) of a 0.5 /im wide plasmon 
transmission line as a function of the frequency for different carrier 
densities at T=77K. 

Fig. 1 101 shows the threshold gain values needed to achieve oscilla- 
tion as a function of the cavity length L for different values of the 
parameter |r|. Comparing the curves in Fig.[TO]to the calculated gain 
values in Fig. |4] it can be seen that the large plasmon gain values in 
graphene could allow plasmon oscillators with cavity lengths much 
smaller than 10 /im. As a practical design example, we consider 
a 0.5 /xm wide, 5 /im long, plasmon transmission line connected 
to an external impedance Z^xt^i^) of 1.0 ¥Sl at 4.0 THz. This 
value of external impedance would correspond to the impedance 
of an integrated dipole antenna plus the impedance of the ohmic 
contacts 1261 . We assume T = 77 K and e^o ~ 4 to. From Fig. |9l 
the impedance Zo(cj) of the transmission line at 4.0 THz is around 
2.2 KQ,. Therefore, r{uj) is -0.63. Using Equation l |30l ), the required 
threshold gain comes out to be approximately 1800 cm~^. This 
threshold gain value is much smaller than the theoretically predicted 
maximum net gain values of around 20,000 cm~^ available at 4.0 
THz in graphene at moderate electron-hole densities (see Fig. ^ 
providing a good margin against extra losses that might not have 
been accounted for in this paper 

The impedance Zs (co) of a graphene strip, given in Equation ( |31t . 
is also interesting since in principle it can be experimentally measured 
at even high frequencies using terahertz time-domain spectroscopy 
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Fig. 11. The real pail of the impedance Zs{w) of a 0.5 wide and 5 
fim long graphene strip is plotted as a function of the frequency for different 
carrier densities at T=77K. The assumed values of v and t are 10* cm/s and 
0.5 ps, respectively. In the presence of net plasmon gain (e.g., for frequencies 
between 2.5 and 6.5 THz for n = p = 3 X 10^" cm^^ - see Fig. |4), the real 
part of the impedance Zs (w) is negative. 
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Fig. 12. The imaginary part of the impedance Zs{lo) of a 0.5 fj,m wide 
and 5 fim long graphene strip is plotted as a function of the frequency for 
different carrier densities at T=77K. The assumed values of v and r are 10* 
cm/s and 0.5 ps, respectively. 



techniques (271. The expression in Equation l l31t indicates that 
the impedance will exhibit resonances in frequency with spacing 
given by the inverse plasmon roundtrip time in the graphene strip. 
These resonances correspond to the confined plasmon modes of the 
graphene strip. Fig. [TT] and Fig. [T2] show the real and the imaginary 
parts of the impedance Zs(io) of a 0.5 /im wide and 5 fim long 
graphene strip as a function of the frequency for different electron and 
hole densities at T=77K. The plasmon resonances in the impedance 
are clearly visible. The resonances disappear when the plasmon loss 
or the plasmon gain are large, and become pronounced when the 
net plasmon gain is close to zero. In the presence of net plasmon 
gain (e.g., for frequencies between 2.5 and 6.5 THz in Fig. [TT] and 
Fig. [12] for n = p — 3 X 10^*^ cm~^ - see Fig. |4}, the real part of 
the impedance Zs(io) is negative. Experimental measurement of the 
impedance Zs (cu) under population inversion conditions can therefore 
be used to measure the plasmon gain in graphene. Note that at very 
small frequencies, when qL >> 1 does not hold, the expression for 
the impedance in Equation ( I31t is not valid. 



6 







/ 

£■=0 = Sj 




y 


X 






m 


graphene -""''^ 

substrate 



Optical pumping 



Fig. 13. Plasmon propagation across a dielectric interface in graphene: (TOP) 
Top view (BOTTOM) Side view. 



VI. Plasmon Confinement and Waveguiding 

In Section ini we mentioned that the width W is not necessarily 
the physical with of the graphene layer. In this Section, we show that 
plasmons can be confined laterally and guided in hybrid waveguide- 
transmission-line structures. This is useful since plasmons need to 
be kept away from regions where population inversion is small in 
order to minimize losses. To illustrate the basic idea we consider a 
graphene layer on top of an insulating substrate. The top surface of 
the graphene layer is covered by two insulating media of different 
dielectric constants such that the value of in the two cases are 
e,; and et, and these two media form an interface at x=0, as shown 
in Fig. [T3] Consider a plasmon wave incident at an angle 9i at the 
interface. The transmitted wave makes an angle 6t. The condition 
for phase matching f25 |, and the continuity of the potential and the 
normal component of the current at the interface can give the value 
of the angle 9t and the amplitudes of the reflected and the transmitted 
plasmon waves. The phase matching condition gives, 



qi{Lj) sm{9i) = qt{ijj) sin(6't) 



(34) 



For small wavevectors, the above Equation can be approximated with 
very little error by. 



ei sin(6'i) = et sm{Ot 



(35) 



If > et, then for 9i larger than a certain critical angle 9c the 
plasmon wave is totally internally reflected. The value of the critical 
angle 9c is sin^^(et/ei). For 9i larger than the certain critical angle, 
the incident plasmon wave is completely reflected. This suggests that 
it is possible to confine and guide plasmons in hybrid waveguide- 
transmission line structures in which vertical confinement is provided 
by the graphene layer itself and lateral confinement is provided by 
the dielectric guide that has a single core sandwiched between two 
claddings on either side. A cross-sectional view of such a guide 
is shown in Fig. [14] The plasmon guide, like dielectric optical 
waveguides |25|, can support more than one transverse mode. A 
single transverse mode plasmon guide is desirable for efficient power 
coupling to the external circuit and can be obtained by keeping the 
difference between the core and cladding too values small. A full 
analysis of the modes supported by the guide shown in Fig. [14] is 
beyond the scope of this paper but we should mention here that for 
a single transverse mode plasmon guide the width W used in the 
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Fig. 14. A cross-sectional view of a hybrid waveguide-transmission line 
structure for plasmon confinement and guiding. Population inversion is 
achieved by optical pumping. 
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Fig. 15. A cross-sectional view of a plasmon guide integrated on a gated 
pn-junction in graphene. Population inversion is achieved via electron-hole 
injection under an applied forward bias. 



analysis is roughly equal to the effective width of the transverse mode 
in the guide. Methods to achieve population inversion in graphene are 
discussed in the next Section. 

VII. Achieving Population Inversion in Graphene 

Schemes for achieving population inversion in graphene via optical 
and electrical pumping have been extensively discussed by Ryzhii et. 
al. in II3I , II4I in the context of obtaining negative conductivity in 
graphene at terahertz frequencies. A knowledge of the electron-hole 
recombination times in graphene is critical for the evaluation of differ- 
ent schemes but no experimental results have so far been reported on 
the recombination times in graphene. Small band-gap semiconductors 
usually have large electron-hole recombination rates due to Coulomb 
scattering (Auger recombination) [28|. The zero bandgap and the 
large optical phonon energy in graphene (> 180 meV |29|) suggest 
that electron-hole recombination rates could be dominated by Auger 
processes. Results in Ref. | 24| show that electron-hole recombination 
times due to Auger scattering in graphene are weakly temperature 
dependent and can range from 1 ps to 100 ps for electron-hole 
densities in the lO^-lO'^^ cm~'^ range. Optical pumping scheme for 
achieving population inversion in a plasmon waveguide is depicted 
in Fig. [14] The electron-hole generation rate Gp in graphene under 
optical pumping with pump intensity Ip and pump frequency up is 
given by the expression |14|, |30|, 



(36) 



, 60 c hujp 

Assuming T=77K, an Auger recombination time of around 10 ps for 
electron-hole densities in the low 10^" cm^"^ range (see Ref. f24^), 
and optical pumping with a 10 wavelength laser, achieving 
plasmon gain at around 4.0 THz would require a pump intensity of 
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2 KW/cm^. This intensity level is easily achievable with a focused 
beam from a standard 5-Watt CO2 laser. High power Q-switched 
pump lasers can also be used at shorter near-IR wavelengths for 
pulsed operation. Another way to achieve population inversion in 
graphene is via electron-hole injection in a forward biased electro- 
statically gated pn-junction, as is done in semiconductor interband 
lasers fill, [131. This scheme is depicted in Fig. [TS] which shows the 
cross-section of a plasmon guide at the junction of electrostatically 
gated p and n regions. The guide helps to confine the plasmons in the 
region where population inversion is expected to be maximum under 
an applied forward bias. The zero-bandgap of graphene suggests that 
band-to-band tunneling current could be large in a forward biased pn- 
junction and could possibly limit electron-hole injection. In Ref. 1 13| 
it was shown that the injected components of the electron and hole 
currents in a forward biased graphene pn-junction could be significant 
and could allow for population inversion to be achieved in the vicinity 
of the junction. 

VIII. Challenges and Conclusion 

In this paper, we proposed and analyzed terahertz oscillator designs 
based on plasmon amplification in graphene via stimulated emission. 
We presented a transmission line model for plasmon propagation and 
amplification in graphene, and obtained the condition for oscillation 
that took into account plasmon losses due to intraband scattering and 
also losses due to external coupling. The ratio of the characteristic 
impedance of the plasmon transmission line to the external impedance 
was shown to play an important role in determining the threshold 
condition for oscillation. The large values of the plasmon gain could 
allow the realization of compact integrated terahertz oscillators. 

Several obstacles need to be overcome before the proposed oscilla- 
tors can be realized. Material quality remains a challenge. Graphene 
monolayers and multilayers produced from currently available ex- 
perimental techniques are estimated to have defect/impurity densities 
anywhere between 10^^ and lO^'^ cm~^ and single-crystal coherence 
lengths of less than a micron 19], f3ll. The values of the electron and 
hole momentum scattering times are critical in determining plasmon 
losses. We modeled the intraband scattering rate with a single time 
constant r the value of which was estimated from experimental 
measurements of carrier mobility in graphene. Several momentum 
scattering mechanisms that could contribute to plasmon losses might 
have not been accounted for in this approach. For example, interband 
momentum relaxing scattering via acoustic phonons could contribute 
to plasmon losses under population inversion conditions for a zero- 
bandgap material. The carrier density dependence of the scattering 
time was also ignored. Plasmon losses due to image currents flowing 
in surrounding metals, ground plane, etc, were not included in the 
analysis presented here. Experimental investigations will be needed 
to provide answers to many of these questions. Given that plasmon 
wavelengths at terahertz frequencies are much shorter than one mi- 
cron (see Fig[3]l, surface roughness of the plasmon guide could scatter 
plasmons and also contribute to plasmon losses. It is also not clear at 
the time of the submission of this manuscript if or not graphene has a 
small bandgap. Several factors, such as disorder, substrate interaction, 
and spin-orbit interaction, can lead to the opening of a small bandgap 
in graphene. The main ideas presented in this paper are expected to 
remain valid qualitatively provided the bandgap is smaller than or 
equal to the frequency of the plasmons. 

The maximum output power levels achievable with the proposed 
oscillators would be determined by the maximum rate at which 
electrons and holes can be injected, device heating, or nonlinear 
effects (such as band-to-band tunneling caused by the plasmon elec- 
tric fields). Nonlinear effects would also become important because 
of the simple fact that the amplitude of the charge density wave 



can never exceed the actual charge density of the electrons and the 
holes in the device. In fact, assumptions made in this paper will 
break down well before the amplitude of the charge density wave 
reaches such levels. A rough estimate of the power levels can be 
made as follows. Assuming an electron-hole density of 6 x 10^*^ 
cm^^, frequency of oscillation of 4 THz, T=77K, a plasmon guide 
of width 0.5 /im and length 5 /im, the values of q{uj) and Zo{i^) 
come out to be approximately 10^ cm^'^ and 2250 Jl, respectively. 
The separation of the electron and hole Fermi levels is approximately 
50 meV. Assuming a density wave amplitude of 2 x 10^" cm~^, the 
amplitude of the potential associated with the wave comes out to be 
approximately 45 mV. This gives a power of the order of 1 /iWatt. 
The compact dimensions of the proposed oscillators could allow small 
phase-locked arrays to increase the output power to mWatt levels. 

The author would like to thank Edwin Kan, Sandip Tiwari, Michael 
Spencer, Faisal R. Ahmad, and Hassan Raza for helpful discussions. 
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